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We calculate the angular dependence of the radiation-zone output power and electric polariza- 
tion of stimulated terahertz amplified radiation (STAR) emitted from a dc voltage applied across 
cylindrical and rectangular stacks of intrinsic Josephson junctions. During coherent emission, a 
spatially uniform ac Josephson current density in the stack acts as a surface electric current density 
antenna source, leading to an harmonic radiation frequency spectrum, as in experiment, but absent 
in all cavity models of cylindrical mesas. We assume that spatial fluctuations of the ac Josephson 
current cause its fundamental mode to lock onto the lowest finite energy cylindrical cavity mode, 
causing it to resonate, leading to a non-uniform magnetic surface current density radiation source, 
and a non-trivial combined fundamental frequency output power with linear polarization for gen- 
eral radiation directions. This combination of the two coherent sources may be fully or partially 
coherent. The coupling of the higher ac Josephson harmonics to other cylindrical cavity modes is 
shown to be very weak, at best, and that emission is also linearly polarized. For the more compli- 
cated rectangular mesas, the lowest energy modes are apparently not excited, but the non-uniform 
ac Josephson current can excite the harmonic sequence of modes with spatial variation across the 
rectangular widths, leading to combined radiation output both for the fundamental and the higher 
harmonics, which combination also may be either fully or partially coherent. We also present a 
model of the superconducting substrate, in which it acts as a perfect magnetic conductor, greatly 
reducing the STAR-emitter power and modifying its angular dependence, especially parallel to the 
substrate. This model can be tested quantitatively by measurements of the angular dependence 
of the emitted radiation, especially that of the second harmonic emitted from cylindrical mesas. 
Quantitative measurements of the amount of coherence could also prove enlightening. Based upon 
this substrate model, existing Bi2Sr2CaCu2 08+(5 crystals atop perfect electric conductors could have 
STAR-emitter power in excess of 5 mW, acceptable for many device applications. 

PACS numbers: 07.57.Hm, 74.50.+r, 85.25.Cp 



I. INTRODUCTION 

The recent discovery of coherent THz radiation emit- 
ted from mesas of the high-temperature superconductor 
Bi2Sr2CaCu208+5 (BSCCO) has caused a great deal of 
excitement [1, 2]. In these experiments, rectangular mesas 
were fashioned by Ar ion milling of a single crystal of 
BSCCO, with a Au layer covering the mesa's top, to 
which an electrical lead was attached, and two additional 
electrical leads were attached to the remaining BSCCO 
crystal substrate. Typical rectangular mesa dimensions 
were 60 x 300 x 1 jiia. By applying a static {dc) voltage 
V across the mesa, the ac Josephson effect was generated 
in each of the N ~ 10^ junctions involved in the mesa, 
and coherent ac Josephson radiation at THz frequencies 
was emitted. The ac Josephson relation 
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where 2e is the magnitude of the Cooper pair charge and 
?i is Planck's constant divided by 27r. Equation (1) relates 
the ac Josephson frequency vj = uoj/{2i:) to the applied 
V/N. To avoid confusion, here we write k and x for the 
wave vectors and positions outside the mesa, and k' and 
x' for the wave vectors and positions inside the mesa. 



Since the frequencies are the same inside and outside the 
mesa, they are left unprimed. 



A. Rectangular mesas 

Since BSCCO is a stack of Josephson junctions, with 
atomically thin superconducting layers separated by 
thicker dielectric layers, it is extremely anisotropic. In 
analogy with standard antenna theory [3-5], to zeroth or- 
der, one might consider a BSCCO mesa to be a dielec- 
tric sandwiched between two metallic layers, forming a 
cavity that produces the electric field that generates the 
radiation [6-20]. In this model, it was usually assumed 
that the ac component of the electric field E{x' ^t) and 
the ac Josephson supercurrent J{x'^t) both had anti- 
symmetric spatial configurations, with maxima along the 
length £ edges, and a node in the center of the mesas of 
width leading to a radiation fundamental wavelength 
A = 2k;, and a finite magnetic field H{x' ^ t) along the ra- 
diating mesa edges, as noted in the following. The wave 
vectors of very thin rectangular cavity modes are 
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where m and p are integers, when there is no spatial vari- 
ation of the electromagnetic fields in the z direction nor- 
mal to the layers. These modes are generated from the 
boundary conditions, which allow for half-integral multi- 
ples of wavelengths across the width and the length of the 
cavity. Usually, such modes in a cavity that is open on 
the sides only arises when the ac magnetic field H{x' ^t) 
vanishes at the edge of the cavity, forcing the normal 
derivative of E{x' ^ t) to vanish there. Other mechanisms 
for the formation of such half-integral modes have also 
been proposed, as discussed in the following. Inside the 
mesa, electromagnetic waves propagate parallel to the 
layers according to the dispersion relation 
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where ^ is the index of refraction, and e ^ 18 for 
BSCCO in the relevant frequency range[l, 2]. By varying 
the dc V experimentally, radiation is found to occur after 
one of the harmonic ac Josephson frequencies nvj locks 
onto a particular (moPo) rectangular cavity mode. 



nUJj = UJmopo 
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In experiments on rectangular samples, it appears that 
the fundamental n = 1 ac Josephson frequency z/j locks 
onto the rectangular cavity (10) mode, so that uj = 
c/{2nrw), also allowing for the higher ac Josephson har- 
monics with n > 1 to lock onto the higher cavity (mO) 
modes, where m = n. 

One might wonder why the low energy excitations were 
not along the length of the rectangular mesa, rather than 
along the width, as observed experimentally. Although 
the experiments on rectangular mesas provided evidence 
for radiation frequencies oc m/w up to four harmon- 
ics of z^j[l, 2], a lower set of modes with harmonic fre- 
quencies would be oc p/£, especially for w <^ £. Very 
recent experiments suggested that standing waves along 
the lengths of rectangular mesas might in fact occur [21], 
although it is presently uncertain if these modes radi- 
ate. A possible resolution of this apparent paradox was 
recently suggested to arise from the magnetic field com- 
ponent normal to the rectangular length, which was pre- 
dicted to have a half- wavelength spatial variation across 
the rectangular width[18]. 

However, this cavity analogy neglects the ac Joseph- 
son supercurrent J{x\t) across each of the junctions. 
Therefore, the second zeroth order model is to treat 
the mesa as a conducting dipole antenna with an ac 
current source, as commonly used in microwave relay 
stations[3, 4, 15, 16, 18, 22]. In this picture, the most 
important part of the ac Josephson current is its portion 
that is uniform across the mesa. This portion leads to 
a non- vanishing H{x\t) within the mesa arising from 
Ampere's law. The important question then, is how to 



treat the Ampere boundary condition. This has been 
done in very different ways, with inconsistent results. 



B. The Ampere boundary condition 

Bulaevskii and Koshelev studied long rectangular 
mesas, and set the magnetic induction B parallel to 
each edge, with E\\z, normal to the substrate [6, 7]. 
They then took \B / E\ = :^(^{kz) at the sample sur- 
face, where (^^ = 1 for lA:;^! < 27rz/j/c = k^, and (^^ 
was imaginary and strongly A:;^ -dependent for kz > k^j. 
The same sort of model for the boundary conditions was 
studied with = Z{uj), a complex constant, by Lin 
and Hu[9-ll], a frequency-independent constant by Lin, 
Hu, and Tachiki[17], and = 7 = 0.1 by Tachiki et 
a/. [18]. Hu and Lin properly took B = on the edge 
of cylindrical mesas, but did not properly take account 
of the ac Josephson current in the Ampere boundary 
condition[13, 14]. Matsumoto et al. treated the Ampere 
boundary condition as providing constraints on B on op- 
posite sides of the long rectangular mesas in their numer- 
ical studies[15, 16]. Koyama et al. treated the ac Joseph- 
son current as uniform in the mesa [20] . The general prob- 
lem with this model for the boundary conditions is that 
a clear understanding of either 7 or Cuj{kz) was not pro- 
vided. In the radiation zone, one expects (^^^ = 7 = ±1, 
but at the edge of the sample, B is given by the Ampere 
boundary condition, which depends upon the ac Joseph- 
son current within the mesa. Thus, the actual situation 
appears at first sight to be much more complicated than 
as suggested in those papers. 

Here we study the simplest and least ambiguous way 
to treat this boundary condition. We use Love's equiva- 
lence principles[3, 4, 23], which are variations on Huy- 
gen's principle, treating the mesa as both an electric 
conductor and as a magnetic conductor [3, 4]. That is, 
the electric field within the mesa is treated as a surface 
magnetic current density source Ms{x\ t) by the electric 
conductor equivalence principle, and the magnetic field 
inside the mesa is treated as a surface electric current 
density source Js{x\t), using the magnetic conductor 
equivalence principle [3, 4]. A sketch of a cylindrical mesa 
with a uniform Js and a nonuniform Ms corresponding 
to the (11) cylindrical cavity mode is shown in Fig. 1. 
This combination of equivalence principles allows us to 
obtain analytic forms for the angle dependence and polar- 
ization of the radiation emitted from cavities suspended 
in vacuum, while properly accounting for the Ampere 
boundary condition. 

As long as rectangular mesas were the only radiation 
sources under study, it has proved difficult to distinguish 
the results of these radiation sources [22, 24]. Clearly, 
the node (or soliton kink) in the center of the mesa as- 
sumed by some to be present in the cavity model cannot 
be the ground state of the system[18, 19], even dynam- 
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Xml Xm2 Xm3 Xm4 Xm5 Xm6 

3.8317 7.0156 10.1735 13.3237 16.4706 19.6159 

1 1.8412 5.3314 8.5363 11.7061 14.8636 18.0155 

2 3.0542 6.7061 9.9687 13.1704 16.3475 19.5129 

3 4.2012 8.0152 11.3459 14.5858 17.7887 20.9725 

4 5.3175 9.2824 12.6809 15.9641 19.1961 22.4013 

5 6.4156 10.5205 13.9872 17.3129 20.5755 23.8036 

6 7.5013 11.7350 15.2682 18.6374 21.9317 25.1840 

7 8.5778 12.9324 16.5294 19.9405 23.2681 26.5450 

8 9.6474 14.1155 17.7740 21.2291 24.5872 27.8893 

9 10.7147 15.2868 19.0046 22.5014 25.8913 29.2186 
10 11.7699 16.4479 20.2212 23.7607 27.1820 30.5345 



FIG. 1: (Color online) Sketch of a cylindrical mesa with a 
surface electric current density Js (red vertical arrows) and a 
magnetic current density Ms (blue horizontal azimuthal ar- 
rows). The =b signs refer to the signs of the (11) cavity mode, 
separated by the yellow line, with fixed angle 00 indicated by 
the thin horizontal black arrow. 



ically, and inhomogeneous, layer-index-dependent kinks 
are even higher in energy [9-11], but the extreme non- 
linear current-voltage characteristic of each mesa allows 
for chaotic, nonequilibrium effects to be important, so 
such complications could potentially play important roles 
in the origin of the radiation. As noted above, addi- 
tional rectangular cavity modes might form across the 
length of the sample [21]. The only practical ways to 
distinguish these models in rectangular mesas is by the 
observations of the angular dependence of the output 
power, polarization, and degree of coherence. As in- 
dicated previously [22], the antisymmetric cavity model 
necessarily leads to a maximum in the radiation output 
power at ^ = 0°, directly above the mesa, whereas the 
conducting dipole antenna model leads to zero radiation 
power at ^ = 0°. The situation near = 90° is much 
less clear. The conducting dipole antenna model leads 
to a substantial output power in that direction, at least 
in the absence of a superconducting substrate. However, 
depending upon the wave vector k of the outgoing ra- 
diation, emission from the cavity model can lead either 
to a vanishing or a non- vanishing output at ^ = 90°, 
at least for a non-superconducting substrate. Thus, we 
considered it important to study a different experimental 
configuration, in which more substantial differences be- 
tween the output predictions of these two models would 
be clearly evident. It is also crucial to consider the effects 
of the superconducting substrate. 



C. Cylindrical mesas 

Here we describe the application of these two mod- 
els to mesas with cylindrical geometry, which recently 
have been studied experimentally [26]. We note that Hu 
and Lin numerically studied the radiation from cylindri- 
cal mesas suspended in vacuum[13, 14, 31]. In place of 



TABLE I: Table of the first six wave vector parameters Xmp — 
kmpCL for m = 0, . . . , 10, where p defines the rank ordering of 



the non- vanishing values of J^iXmp) 
cavity of radius a. 



0, for a cylindrical 



Eq. (2), the mode frequencies for a cylindrical cavity of 
radius a with the boundary condition that the azimuthal 
component of the magnetic field vanishes at the edge, are 
given by 
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where the anharmonic Xmp = ^'mp^ values are listed in 
Table 1. From fits to the fundamental mode on three 
cylindrical mesas [26], the choice rir ~ a/TS fit the as- 
sumption that the fundamental ac Josephson frequency 
locked onto the cavity (11) mode reasonably well, as it did 
in fitting the rectangular cavity (10) mode[l, 2]. Hence, 
in our numerical calculations, we shall use rir = vTs for 
the frequencies of interest. 

Although a layered superconductor consists of a stack 
of Josephson junctions that in some circumstances may 
behave rather independently [25], in our model, we as- 
sume that under the application of a dc voltage V across 
the N layers, all of the layers behave identically. This as- 
sumption was also made by Tachiki et a/., [18], and is con- 
sistent with the conclusions of Bulaevskii and Koshelev 
[6-8], but is distinctly different from that made by Hu 
and Lin[9-14]. For the electric dipole antenna model, 
the frequencies of the emitted radiation are the same as 
those present in the ac Josephson current, Un = nuj^ 
since the integer harmonics are generated by the nonlin- 
ear ac Josephson effect. For the cavity model, however, 
the frequencies Ump of the radiation are the non- vanishing 
values generated by the p zeroes of the first derivative of 
the mth regular Bessel function Jm{z), which are far from 
integer multiples of one another, as evidenced in Table 1. 
If a particular cavity mode were to become in resonance 
with the fundamental ac Josephson frequency i/j, then 
higher harmonics of the Josephson current would not be 
resonant with any other cylindrical cavity mode. Hypo- 
thetically, if some cavity modes did not require resonance 
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with the Josephson current to radiate, the frequency 
spectrum of the emitted radiation would be distinctly 
nonharmonic, as pictured by Hu and Lin [14]. Hence, if 
resonance of the fundamental ac Josephson current fre- 
quency with a non- uniform cavity mode were crucial for 
the radiation to occur, one would not expect to observe 
any higher harmonics in the output power. Observation 
of higher harmonics would therefore be prima facie ev- 
idence that a substantial amount of the radiation does 
not arise from the excitation of cavity modes. 

Here we calculate the angle and frequency dependence 
of the output power for both the cavity and electric dipole 
models for cylindrical and rectangular mesas. We show 
that the experiments on cylindrical mesas presented by 
Kadowaki et al. provide evidence that the radiation at 
the fundamental ac Josephson frequency, consistent with 
the (11) cylindrical cavity mode, is a mixture obtained 
from both the electric and magnetic current sources [26]. 
To fit the data quantitatively, we must also take account 
of the dramatic effects of the superconducting substrate. 
However, the higher frequencies observed in the measured 
output power are harmonics of the fundamental, so that 
any mixing with higher cavity modes is very unlikely. 
Measurements of the angular dependence of the output 
power, polarization, and coherence of the second har- 
monic of the emission from cylindrical mesas atop super- 
conducting substrates should therefore be fully coherent, 
have linear polarization, and vanish at ^ = 0° and 90°. 
Studies of the angular dependence of the power, polariza- 
tion, and coherence of the emission at the fundamental 
frequency may help to separate the two radiation sources, 
as well. For rectangular mesas, the correct boundary con- 
dition oi B\\ = is difficult to impose precisely in a rect- 
angular geometry, due to the corners. However, with the 
Love equivalence principles, one can find a closed form 
expression for the far-field radiation, the results of which 
will be presented, as well. 

The paper is organized as follows. In Sees. H-VHI, the 
formation of the radiation emitted from cylindrical cavi- 
ties is studied theoretically, since the mathematics is the 
simplest. In Sec. II, we focus on the primary radiation 
source, which we take to be the ac Josephson current. 
In Sec. Ill, we study the spatial fluctuations within the 
mesa of this primary radiation source. In Sec. IV, we de- 
scribe how these spatial fluctuations may excite a cavity 
mode. In Sec. V, we discuss the radiation from a cavity 
mode alone. In Sec. VI, we calculate the output power of 
the combined primary and excited mode secondary radi- 
ation. In Sec. VII, we calculate and discuss the electric 
polarization and coherence of this combined radiation. In 
Sect. VIII, we discuss and model the effects of supercon- 
ducting substrates. In Sec. IX, we present our analogous 
results for rectangular mesas. We discuss and summarize 
our results in Sec. X. 



II. PRIMARY RADIATION SOURCE 

As for rectangular mesas, a dc voltage V is applied 
across the stack of N superconducting junctions in the 
BSCCO mesa, which causes the junctions to exhibit an ac 
Josephson current, which oscillates at the frequency vj 
given by Eq. (1). In addition, the intrinsic non-linearity 
of the Josephson junctions causes the ac Josephson cur- 
rent to have a large number of harmonics at Vn = nvj^ 
regardless of the spatial dependence of the Josephson cur- 
rent within each junction. From Love's magnetic equiv- 
alence principle, the magnetic field in the cavity gener- 
ated by the ac Josephson current along z according to 
Ampere's law is treated by replacing it with a surface 
electric current J^, and setting the resulting magnetic 
field in the cavity equal to zero [3, 4]. This ac Josephson 
current has two essential functions in the radiation. First, 
it radiates at all of its harmonic frequencies. Second, the 
radiation at one of its frequencies may lock onto a cylin- 
drical cavity mode, exciting it, and causing it to radiate, 
as well. Hence, the ac Josephson current is the primary 
radiation source. We assume that during emission, all of 
the junctions radiate together, so we neglect the layer 
index of the spatial variation [2 5], and write within the 
mesa 

oo 

J{x',t) = i'^e-»--*[J„^ + <5J„(a;')], (7) 

n=l 

where the spatial average {8Jn{x')) in the mesa vanishes. 
We assume that any time dependence to 5Jn{x') is slow 
with respect to the measurement times, and can be ne- 
glected. 

We then may calculate the resulting radiation pattern 
in the conducting electric dipole antenna model by taking 
the equivalent electric current source to be the integrated 
current source placed on the surface from which the ra- 
diation occurs, which is the edge of the cylindrical mesa 
sketched in Fig. 1. We assume no radiation emanates 
from the top and bottom surfaces of the cylindrical mesa. 
This leads to the surface electric current density, which 
we write just inside the mesa edge as 

Js{x',t) = ^rj{zW-^)J{^'.t), (8) 

where J{x\t) is given by Eq. {7),r]{z') = Q{z')Q{h- z') 
is unity over the height h of the mesa, a is the radius of 
the cylindrical mesa, and we assume /i/a <C 1. In this 
limit, it suffices to take r]{z') h5{z') when the mesa is 
suspended in vacuum. A sketch of the uniform part of 
the surface electric current density Js at the edge of a 
cylindrical mesa is given in Fig. 1. 

If the electric surface current density contained only 
the fundamental ac frequency i/j, the magnetic vector 
potential Ai(cc, t) arising from it would be given by 
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where R= \x — x'\. We note that uoj is the same inside 
and outside the mesa, but kj = ujj/c outside the mesa 
must match the wave vector of hght in vacuum. Inside 
the mesa, kj — n^cjj/c. Since A\\z is parahel to the 
mesa edge, this wave vector automat icahy satisfies the 
Maxwell boundary condition that the tangential compo- 
nents of E are preserved across the boundary. General- 
izing this standard expression to the full set of harmonic 
ac Josephson frequencies, we have 



X[J^ + 5Jn{x')]. 
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We write 



A{x,t) = Aj{x,t) + Asjix,t), 



(10) 



(11) 



separating it into its contributions from the spatially uni- 
form and inhomogeneous 5J:^{x^) source amplitudes. 
In Sec. Ill, we present a theory for 5J{x')^ and use it to 
evaluate A5j{x^ t). 

In the far field, or radiation zone, r /a ^ 1, and we use 
the standard approximation 
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where \x\ = r in spherical coordinates [3-5]. In the radia- 
tion zone, z' —6 sin 6 and k x' = kp' sin cos(^ — (j)'). 
Taking h/a <C 1, so that ri{z') h5{z') when the mesa 
is suspended in vacuum (in a Gedanken experiment), the 
integral in Eq. (10) is then readily found to be 
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where v = ira'^h is the volume of the mesa, kg = kja sin ^, 
and S:^{0) = 1 in the absence of a substrate. Then, it is 
elementary to obtain EA{x^t) from Maxwell's equation, 
Ea = since there is no electric potential outside 

the mesa. Similarly, Ha = x A. In the radiation 
zone, the contribution Eaj from Aj is given by 



EAj{x,t) 



iO sin OvfiQ 
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xnJiujjJ^{nke)Si{e). (14) 
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FIG. 2: (Color online) Plots of the radiation intensity in arbi- 
trary units emitted from the uniform part of the ac Joseph- 
son current when a cylindrical mesa is suspended in vacuum, 
(a) The radiation at the fundamental n — 1 mode, (b) The 
radiation at the n — 2 second harmonic mode. 



We note that this implies the radiation from the sur- 
face electric current source is linearly polarized, with its 
polarization direction parallel to ^. So far, however, 
kj = cjj/c is only determined by the applied V/N ^ and is 
unspecified with respect to the mesa dimensions. When 
the fundamental ac Josephson mode with frequency tjj 
locks onto the imp) cavity mode, the resulting k'j and 
kj values are given by k'ja = uojUr/c = Xmp inside and 
kja = uoj/c = Xmp/^r outside the mesa, in order to 
properly describe the electromagnetic propagation in the 
vacuum. 

The radiation-zone intensity /(^, (j)) for the conducting 
dipole model is given by the differential power per unit 
solid angle by dP/dQ = ^Re[r^f • E x i^*][5], where — 
is a time average. The contribution from the uniform 
source in the radiation zone is 



^ r/7»i Sin' 0Y,\Br.{0)Mnke)\^ (15) 



w here B njO) = J^vS^{0)nkj^/Z^/{4V27r), and Zq = 
Y^/io/^o is the vacuum impedance. The radiation pat- 
terns expected with kja = 1.8412/nr = Xii/^r from the 
uniform part of the conducting dipole model alone at 
the fundamental (n = 1) and second harmonic (n = 2) 
when the mesa is suspended in vacuum are pictured in 
Figs. 2(a) and 2(b), respectively. These intensity pat- 
terns are nearly indistinguisable, except for the emission 
frequency, due to the long wavelength in vacuum rela- 
tive to the mesa size. The intensities vanish at ^ = 0°, 
due to the geometric factor sin^ in the radiation-zone 
spherical coordinate representation of z', are indepen- 
dent of (/), and are non- vanishing at ^ = 90°. When a 
cylindrical mesa is suspended in vacuum, /(^, cj)) for the 
fundamental n = 1 and second harmonic n = 2 modes 
each exhibit a maximum at ^ = 90°, as for the ordinary 
dipole model based upon a quasi-one-dimensional wire or 
thin rod-shaped source. 
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III. SPATIAL FLUCTUATIONS OF THE ac 
JOSEPHSON CURRENT 

If the ac Josephson current were independent of posi- 
tion within the mesa, this could only couple to a uni- 
form electric field, which would not radiate. Hence, 
the non-uniform part 5J{x') of the ac Josephson cur- 
rent must provide the coupling. There may be a number 
of sources for this non- uniformity, such as defects and 
thermal fluctuations, which may be aggravated by inho- 
mogeneous heating effects [21]. Here we assume that the 
primary source of the inhomogeneities is thermal fluctu- 
ations, but our model can also be employed to approx- 
imate smooth stoichiometry variations, as well[27, 28]. 
In the samples used in the experiments [26], the BSCCO 
crystals were prepared slowly over a very limited central 
region in gold-plated elliptical heating ovens, with highly 
controlled heating, producing samples that are proba- 
bly much more uniform chemically (although still highly 
non-stoichiometric) than in those studies[29, 30]. 

Assuming the magnitude of the fluctuations in the 
spatial dependence of the ac Josephson current is small 
within the mesa with respect to its average value, the 
free energy of these fluctuations may be written as 



Tf oc 



)' + (C)'(V'W„)2], (16) 



where ^^(T) is a temperature-dependent characteristic 
length over which the spatial fluctuations in ^J^ occur, 
and n describes the intralayer spatial variations asso- 
ciated with the nth harmonic of the ac Josephson fre- 
quency, as noted in Eq. (7). Tf is a minimum when 



(17) 



which is written in cylindrical coordinates. The full solu- 
tion to this equation is elementary, and may be written 



as 



m=0 



(18) 



where the 0on and Cm^ are free parameters. Of course, 
we require that the ac Josephson current vanishes out- 
side the mesa (p' > a), but since the uniform part 
is discontinuous at the mesa edge, SJn may also be dis- 
continuous dit p' = a. Instead, we assume that thermal 
fluctuations are weak near the mesa edge, so that we take 

dp' p'=a 

With this von Neumann boundary condition, SJn{x') 
may be written as 



0. 



(19) 
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SJn{x) = C^^Jm{k^pP) 
p=l,m=0 



cos[m(0' - 0on)], 



(20) 



where Xmp = ^mp<^ is the pth non- vanishing value of 
dJm{x)/dx = 0, and we write k'^^ = 1/^^p. The Xmp 
values for m = 0, . . . , 10 and p = 1, . . . , 6 are give in Ta- 
ble I. Since xdxJo{kopx) = {a/kQ^)Ji{kQ^a) = for 
all allowed p values, the spatial average of 5Jn{x') van- 
ishes. We note that k'^^ represents the set of 1/^^ values 
that satisfy our assumed boundary condition. In addi- 
tion, since and SJnix') are real, all of the Cmp are 
real. 

In the Appendix, we calculate the electric field 
EAsji^^t) due to this inhomogeneous ac Josephson cur- 
rent distribution and the resulting combined EA{x,t) = 
EAj{x^t) -h EAsj{x,t), and discuss the complications 
that arise. 



IV. EXCITATION OF A CAVITY MODE 

We next consider the cavity model of the short cylin- 
drical mesa. Since a voltage is applied across the height 
of the mesa, the magnetic vector potential Az satisfies 
^''^Az — p^e^Q^ = —pJz inside the mesa, where p and 
e are the magnetic permeability and dielectric constant 
inside it, and rir = ^/pe ~ y/e. In the absence of an 
ac Josephson current, = 0, the homogeneous solution 



IS 



m=0,p=l 

X cos[m(^' - ^o)], (21) 



where (j)^ and the Bmp 



are arbitrary constants, and we 
neglected spatial variations along the z axis. As part of 
Love's magnetic equivalence principle, Eq. (21) is subject 

to the condition = = [3, 4]. The 

^ ^ p' —a p'—a 

consequence of this requirement and the related feature 
of setting H{x\t) = inside the cavity, is to replace 
the internal magnetic field by a surface electric current 
density Js- Hence, this boundary condition is precise, 
and is not an assumption, as we made for SJn- This 
condition then leads to the same wave vectors k'^^ = 
Xmp /a as for SJn- Adding the spatially inhomogeneous 
part SJn{x') to the wave equation for Az^ we have 



V'^A 



pe- 



—inujjt/^{n) 
c ^mp 



m=0;p,n=l 
y^Jm{k' p') COs[m(0' - 0On)], 



(22) 

which can be solved with the particular resonance solu- 
tion form 

OO 

A^^^ {x ^t) = ^ ^ 6 \.^nipif)^'rn,moSp^pQ 

m=0;n,p=l 
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provided that 

.2 



iC 



mp 

X cos[m 



,(n) 



</>0n)], 



(23) 



2 2 



"mp 



where we have set k'rnopo 



for (mp) = (moPo),(24) 
for (mp) = (mopo), (25) 

for (mp) 7^ (moPo),(26) 

loPo^r/c. The spatiahy 
constant part of Jj{x' ^ t) cannot excite any cavity modes. 
Although Eqs. (24) can be solved for general (mp), Eq. 
(25) can only hold for one value of (mp), which we denote 
(mopo). Of course, if no (mopo) satisfies Eq. (25), then 
none of the cavity modes will be in resonance with any 
of the modes of the inhomogeneous ac Josephson cur- 
rent. Hence, we conclude that the non-linear, spatially 
inhomogeneous Josephson current can excite at most one 
electromagnetic cavity mode, and the other modes are 
off resonance, with Lorentzian lineshapes as suggested 
by Eq. (26). Of course, effects such as a built-in z- 
dependence to the cylinder radius a during mesa fab- 
rication will give finite widths ex 1/Qmp to the cavity 
modes. Note that the t in the solution for the excited 
{mopo) mode causes the amplitude A^lpo {t) of that mode 
to grow linearly in time, which continues until satura- 
tion is achieved, and radiation occurs. During radiation, 
the amplitude of the resonant frequency is approximately 
constant in time, which we take to have the finite value 



l(n) 



,(oo) = lim^ 



l(n) 



,(t). We note from Eq. (24), 

l(n) 



that for small times, the resonant amplitude Amopo {t) is 
pure imaginary for a real e in the frequency range of in- 
terest, so we assume that it stays pure imaginary after 
saturation, and write 



^mopo 



(oo) 



^mopo 

2neujj ' 



(27) 



where teff is the effective time to reach saturation. 

Since the lowest energy cavity mode is the (mp) = (11) 
mode, if all of the junctions are involved in the radiation, 
it will be the easiest mode to excite by experimentally 
increasing V. If the upturn in the current-voltage (I/V) 
characteristic occurs before this mode can be excited, 
then no resonance of the ac Josephson fundamental mode 
including all N junctions with any cavity mode can be 
achieved. 



V. RADIATION FROM A CYLINDRICAL 
CAVITY MODE 

Now, we focus our attention on the angular depen- 
dence of the radiation, assuming that a particular cavity 



n=1 



n=2 n=3 n=4 



n=5 



(11) 



(21) 



(01) (41) 
)(31)(12) 



(02) 
(22)1 
(51) 



(32) 
(71) 
(611) (13) 



(91) 
(42) (52)/ 
(23) J / 
81X03)/ 







6 
ka 



8 



10 



FIG. 3: (Color online) Sketch of the lowest five wave vectors 
of the ac Josephson current harmonics (long thin black) and 
the lowest 19 cylindrical cavity modes (short thick red), as- 
suming the fundamental (n = 1) ac Josephson frequency locks 
onto the lowest energy (11) cavity mode frequency. The (mn) 
indices of the cavity modes are indicated. 



mode {mp) = (moPo) is involved in the radiation. We as- 
sume that after equilibrium has been achieved, the mag- 
netic vector potential from which the electric field in the 
cavity can be calculated is given by the particular form, 
Eq. (23). Then, the magnetic surface current density 
Ms{x'^t) = —n X E^P\x'^t), evaluated at the sample 
edge, [3, 4] where h = p' ^ 



Ms{x',t) 



iip) 



{x',t) 



(28) 



where r]{z') is defined following Eq. (8). 

For a single frequency co and wave vector /c, one obtains 
the electric vector potential F{x,t) outside the sample by 
the standard integral 

F{x,t) = ^ I (fx'Ms{x',t)e'^^/R, (29) 
47r J 

where Ms{x'^t) is the surface magnetic current density 
at the frequency co = c/c[3-5]. 

The generalized F{x^ t) is then obtained from Eq. (29) 
by replacing uo with oo^nopo = ^^j, and the details are 
given in the Appendix. Using the Schelkunoff expression 



for the electric field, Ep = — — V x F, we obtain for the 



_ _x 

resonant contribution to Ep in the radiation zone 

) ^ 

x{0cos6>sin[m((/) - 0On)]-^m(^mp) 



E 



F (^7^) r/a>l ( ^) ^rnl' 



8 



0=0 




FIG. 4: (color online) Plots of the radiation intensity in ar- 
bitrary units emitted from cavity modes when the mesa is 
suspended in vacuum, (a) The cavity (11) mode with kna = 
1.8412/nr. (b) The cavity (12) mode with ki2a = 5.3314/nr. 



0-0 



(a) 




FIG. 5: (Color online) Plots of the radiation intensity in ar- 
bitrary units emitted from cavity modes when the mesa is 
suspended in vacuum, (a) The cavity (01) mode with koia = 
3.8317/nr. (b) The cavity (21) mode with k2ia = 3.0542/nr. 



The magnetic field Hf at the frequency of interest may 
be obtained from = V x Ep. We note that the 

overall factor — (— i)^ in E^^^ plays a very important role 
when Ea is added to as discussed in the following. 

We then calculate the output power of the (moPo) cav- 
ity mode in the radiation zone, which is averaged over t 
after the resonance has saturated. We note that some or 
all of the (j)on might vary with experimental run, in which 
case an average of the output power over all possible (j)on 
values might be warranted. Assuming the (moPo) mode 
of the cavity model of a cylindrical mesa can be excited 
sufficiently in order to radiate, its differential power per 
unit solid angle in the radiation zone at fixed (j)on is then 



mp 



327r2 



-(^cosl9sin[m(0 - (l)On)]J:^{ktp) 



2i 



(31) 



In the radiation zone, spherical plots of the radiation 
patterns for the four TM^^ cavity modes with the lowest 
energies and with 0on fixed are shown for a cylindrical 



mesa suspended in vacuum in Figs. 4-5. For fixed ^on, 
the radiation patterns all exhibit C2mv point group sym- 
metry, invariant under rotations of 7r/m about the ^ = 
axis, and exhibiting 2m mirror planes. For m = 0, we 
may denote the point group symmetry as Coov^ as it is 
rotationally invariant about ^ = 0°, with a limitless num- 
ber of mirror planes. We note that the far-field radiation 
from the m = 1 cavity modes, two of which are pictured 
in Fig. 4, exhibit absolute maxima at 6> = 0, suggestive of 
a relationship with the experiments, and they exhibit a 
substantial azimuthal anisotropy, especially for ^ 90°, 
where the output from the experiments is small or may 
even vanish [26], unlike the azimuthal average of the pat- 
terns shown in Figs. 4 and 5. The output power when the 
mesa is suspended in vacuum from all other cavity modes, 
such as the (01) and (21) modes pictured in Fig. 5 and 
those of the not-pictured (31) and (41) modes, vanish at 
^ = 0, unlike the experiments. All cavity modes emitted 
from mesas suspended in vacuum have large outputs for 
6 = 90°. The output power of the cavity (01) mode pic- 
tured in Fig. 5(a) is independent of ^, as is evident from 
Eq. (31), and is very similar to that of the ac Josephson 
n = 2 mode pictured in Fig. 2(b). We remark that Hu 
and Lin presented very similar figures for the same four 
modes that we showed in Figs. 4-5 recently [14, 31]. 



VI. COMBINED PRIMARY AND SECONDARY 
RADIATION 

After equilibrium has been achieved, we assume the 
fundamental n = 1 ac Josephson mode has locked onto 
the cavity (moPo) = (H) mode, with the wave vector 
kiia = 1.8412/nr, as seen in experiments on three cylin- 
drical mesas [26]. The combined electric field will then 
have the components 



Ee 



i[G^^^ cos((/) - (/)oi) Ji [ke] - Bi sin6>Jo(/c^)] 

^i{kjr-ujjt) 



iC'ii cosOsm{(j) — Ji {kg) - 
nviiQUJjJiSi{e), 



i{kjr-ujjt) 

47rr 



(32) 

-,(33) 
(34) 



Because — (— i)^ = z for m = 1, the components to Eq 
from the spatially homogeneous J:^ contribution to Ea 
and from E'p^ are either in phase or out of phase by tt, 
depending upon ^oi- However, when the inhomogeneous 
contributions 5J^{x') are included in Ea^ as shown in 
Eq. (50), the contribution proportional to cos[m(0 — 
</>0n)] to the full Ea is proportional to (— z)^+^, so that 
this contribution would be out of phase by ±7r/2 with 
the two terms included in Eq. (32). If significant in 
size, such terms would have an important, experimentally 
observable effect, both upon the angular dependence and 
the polarization of the radiation. 



9 



The combined output power of the fundamental n = 
1, spatiahy homogeneous ac Josephson surface current 
density mode locked onto the cavity (11) mode at coj = 
coii is given in the radiation zone by 



dP^ 



11 



^ r(G^Vcos^sin((/)-(/)oi)J+(fe^)) 
BismOJo{ke)y\. (35) 



Note that since the contributions to Eq from electric and 
magnetic surface current densities are either in phase or 
TT out of phase for the (11) cavity mode locked onto the 
n = 1 dipole mode, the fully coherent output power is 
asymmetric in (/), breaking the C2v point group symmetry 
to Ci. There remains only a single mirror plane, which is 
normal to the substrate and contains ^oi- A plot of the 
combined output power from the (11) and n = 1 modes 
for Bi/G'ii = 5 when ^oi is fixed and the mesa is sus- 
pended in vacuum is shown in Fig. 6(a). This figure is in 
rather good agreement with experiment, except for the 
behavior near to ^ = 90°, where the experimental values 
are very small, and may even vanish [24]. There is a point 
kink, a point where the output power is much smaller 
than in the surrounding solid angle, which lies along the 
^01 line. Detector scans across the cylindrical mesa top 
that pass this point are only infinitesimally probable, and 
hence may record a significant output power when the 
detector passes nearby but not directly over the position 
of the point kink. However, we note that this simple 
output form arose from the neglect of the spatially inho- 
mogeneous ac Josephson current SJ^{s') source to Ea, 
as discussed in Sec. III. Inclusion of the contribution 
proportional to cos{(j) — (j)oi) would be =b7r/2 out of phase 
with the two components in Eq. (32), and would tend 
to remove the point kink present in Eq. (35) and shown 
in Fig. 6(a), and weaken the overall degree of coher- 
ence. Higher azimuthally anisotropic terms proportional 
to cos[m(^ — ^oi)] would also be present. Hence, if it can 
be established both that ^oi is not random, and second 
that such a kink in the angle dependence of the out- 
put power is absent, then that would provide evidence 
for a second role for the SJ^{x') inhomogeneous electric 
current source, besides its excitation of the appropriate 
cavity modes. Note that the output power is finite at 
both i9 = 0° and i9 = 90°. 

In Fig. 6(b), a plot of the output power of the incoher- 
ently combined n = 1 dipole and (11) cavity modes with 
001 random and Bi/G^i = 7 is shown. This output is 
taken by averaging Eq. (31) over (/)oi, which can be done 
by inspection. But it has a drastic effect on the output 
point group symmetry, which is now Coov , exhibiting ro- 
tational invariance about ^ = 0°. As seen in Sec. VIII, 
it also has a strong effect on the polarization. 

From the data in Table I and the graph in Fig. 3, it 




FIG. 6: (Color online) Plots of the radiation power emitted 
from the mixed n = 1 ac Josephson and (11) cavity (11) mode 
with 001 fixed and kna = 1.8412/nr for the mesa suspended 
in vacuum, (a) Bi/G^^^ = 5 and is ^oi fixed as indicated, (b) 



Bi/G^ 



(1) 



7 and ^oi is random. 



is evident that when the n = 1 fundamental ac Joseph- 
son mode locks onto the cavity (11) mode {kj = k[i = 
kii/rir)^ the cavity (01) mode wave vector k^^ only ex- 
ceeds that of the second ac Josephson harmonic 2kj by 
^ 4 %. Although it will be nearly impossible to excite 
both the cavity (11) and (01) modes simultaneously, in 
case the Qoi of the cavity (01) is very low, that mode 
might be very slightly excited, slightly affecting the out- 
put power of the radiation at 2a; j. The combined electric 
fields at 200 J = 2k ja after equilibrium has been achieved 
are given by 



Ee 



-iB2 smOJo{2ke) 

i{kjr-ujjt) 



(36) 



where B2 and Gqi are given by Eqs. (34) and (56) in the 
Appendix. Note that the off-resonant cavity response to 
the driving ac Josephson nth harmonic differs in phase 
by — 7r/2 from the resonant response, as noted just before 
Eq. (56). Hence, the contributions from the second ac 
Josephson harmonic current source and the off-resonant 
(01) cavity mode are out of phase by 7r/2. We note that 
the combined radiation exhibits Goov point group sym- 
metry, rotationally invariant about ^ = 0°, regardless of 
whether (j)oo is fixed or random, and is linearly polar- 
ized along 0. Since both terms vanish at ^ = 0° and we 
expect |Gq^7^2| <C 1, the output power is unlikely to 
differ noticeably from that pictured in Fig. 2(b), except 
for the substrate factor discussed in Sec. VIII. We note 
however, that some azimuthal anisotropy to the second 
harmonic output power could arise from the spatially in- 
homogeneous part 8J:^{x') of the ac Josephson current, 
as shown in Sec. HI. From Eq. (50), the n = 2 contribu- 
tion could have various non- vanishing terms proportional 
to cos[m(0 — 001 )] 7 although we expect their amplitudes 
to be small. 
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VII. ELECTRIC POLARIZATION OF THE 
COMBINED RADIATION 

The polarization of the combined output at the funda- 
mental frequency uoj = cJn and kj = /cn/n^, under the 
assumption that the conditions are appropriate for the 
fundamental Josephson frequency to excite the (11) cav- 
ity mode is more interesting. After equilibrium has been 
achieved, the combined electric field in the radiation zone 
will have the components given in Eqs. (32) and (33). 
The radiation of the combined Josephson n = 1 funda- 
mental mode locked onto the cavity (11) mode generally 
has linear polarization. This is because the amplitudes 
of Eq and E(fy are pure imaginary, so the phase differ- 
ence between these two components is zero for all times. 
However, the axis of polarization depends upon the par- 
ticular values of (6>,0). There are some special cases. 
These are (I) If is random on the time of the measure- 
ment, the radiation is linearly polarized along ^, arising 
from the dipole radiation only. This vanishes at ^ = 0°, 
of course, so the radiation would be unpolarized there, 
and one would measure no radiation power through the 
polarizer. (II) If ^ — ^oi = 0, tt, the polarization is along 
0. (Ill) If ^ = 0°, the polarization direction depends 
strongly upon (j) — But in a measurement at 6> = 0°, 
one would integrate over all (j) values, measuring an un- 
polarized radiation, in agreement with experiment [24]. 
(IV) If 6> = 7r/2, the radiation is linearly polarized along 
0. (V) At the special line where Eq = {)^ the polarization 
is along 0. There is in addition one line given by 

where Eq = 0, and the polarization is along cj). 

For the general case, the polarization is tilted away 
from the 6 direction by the angle r given by [3, 4] 




In Fig. 7, we presented three-dimensional plots of 
2r/7r in radians versus 6 in degrees and r = Bi/G^i for 
^— (/>oi = 7t/2 and (j)—(j)oi = 7r/4 in the left and right pan- 
els, respectively. In Fig. 7(a), when cj) — 0oOl = 7r/2, Eq 
is only given by the ac Josephson (or dipole) radiation, 
which vanishes in the limits ^ = 0° and r = 0. When 
= 90°, as long as r > 0, the polarization is along 0. 
From this figure, it is evident that when r ~ 1 or greater 
and 6 ^ 10° or greater, r 0, and the polarization lies 
along 6. In Fig. 7(b), when ^ = 0° or r = 0, the polar- 
ization is -45° between cj) and 0^ which is clear from Eqs. 
(32) and (33). Then, as either r ~ 1 and ~ 10°, the 
polarization switches to the 6 direction, assuming the co- 
herent combination with 0oi is fixed in the measurement. 




FIG. 7: (Color online) Plots of the normalized tilt angle 2r/7r 
in radians that the polarization makes with 6 versus 6 in 
degrees and r — Bi/G^^^ of the mixed source radiation for 
fixed 001 at the fundamental n — 1 frequency with kna = 
1.8412/nr when the mesa is either suspended in vacuum or 
on a superconducting substrate, (a) — 0oi = 7r/2. (b) 
— 001 — 7r/4. . 




FIG. 8: (a) Sketch of a mesa with Js and Ms surface current 
sources, (b) Mesa during coherent emission with applied dc 
/o, Vb, and lac confined to it. Curve c is the integration path 
for the Ampere law boundary condition. See text. 



If 001 is not fixed, but random, then the polarization will 
always be along ^, except that one will measure an unpo- 
larized radiation at ^ = 0°. Preliminary measurements 
confirmed the unpolarized radiation at ^ = 0°. However, 
that single experiment cannot distinguish whether 0oi is 
fixed or random, because at 6> = 0°, the polarization av- 
erages to zero, both for fixed and random 0oi. Hence, 
further experiments are necessary to distinguish whether 
001 is fixed or random. A measurement of the degree of 
coherence could prove useful in this regard. 



Electric Magnetic Sources 

I — ^ I ^ actual 

I — ► A ^ PMC 

Y J image 

FIG. 9: (Color online) Sketches of the electric (dashed, black) 
and magnetic (solid, red) actual and image sources and on a 
perfect magnetic conductor substrate. From Ref. (2). 
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VIII. SUPERCONDUCTING SUBSTRATES 

We now discuss the consequences of mixed sources of 
the radiation. If a cavity mode were the only source of 
the radiation, the output power would vanish at 6> = 90°. 
This is in close agreement with experiments on rectangu- 
lar and cylindrical mesas [24, 26]. However, if the funda- 
mental ac Josephson current is also a radiation source, it 
would have a contribution at ^ = 90° that would be non- 
vanishing. Since our qualitative fits to the data on both 
types of mesas suggest that the majority of the radiation 
source arises from the electric surface current (or the ac 
Josephson current), it then is important to discuss why 
the experiments see a maximum output at ^ ~ 30°, but 
a seemingly vanishing output power at ^ = 90°. 

Here we argue that the ac Meissner effect in the super- 
conducting substrate sketched in Fig. 8 (a) causes it to 
behave as a perfect magnetic conductor (PMC) [3, 4, 22], 
which can be treated by the image source technique as 
sketched in Fig. 9. Note that in our treatment of the Love 
equivalence principle, we assumed that the magnetic con- 
ductor model allowed us to represent the magnetic field 
internal to the mesa as a surface electric current, and to 
set H = on the surface of the mesa. Present mesa 
antennas atop a PMC substrate emit dramatically lower 
output power than those atop a substrate that is an in- 
sulator or perfect electric conductor (PEC) [3, 4]. We 
note that a superconducting Nb substrate was previously 
used as a ground plane for a square lattice of Nb/Al2 03 
Josephson junctions, the output of which was amplified 
by a cavity at the edge of the array [32-34]. Without the 
ground plane, the radiation was not detectable. In that 
case, the ac Josephson current source was parallel to the 
substrate [32-34], so that the image of the supercurrent 
in the magnetic conducting substrate was parallel to the 
source current, as sketched by the two arrows second from 
the left in Fig. 9. In our case, both the electric and mag- 
netic image currents are opposite to the respective source 
currents, as in the far left and far right arrow pairs in Fig. 
9, causing a cancellation of the output power at 6> = 90°. 
For pure cavity mode radiation at the precise mode fre- 
quencies, the output would vanish there anyway, but for 
the electric current {ac Josephson) source, it would not. 
A study of the angular dependence of the second har- 
monic would therefore be a critical test of this substrate 
model. 

This unexpected combination of 7(90°) ^ and 
^max ^ 30° led us to consider the substrates of existing 
samples [1, 2]. Each existing mesa is a structure formed 
by ion milling on the top center of a wider and thicker 
single crystal of superconducting BSCC0[1]. We neglect 
the thin top Au mesa contact. As sketched in Fig. 8 (b), 
during coherent Josephson radiation, the ac Josephson 
current is essentially confined to the mesa by the applied 
dc current Iq and voltage Vq. With only a dc surface 




FIG. 10: (Color online) Plots of the radiation power emitted 
from the mixed n = 1 ac Josephson and the (11) cylindrical 
cavity modes with kna — 1.8412/nr fixed and the mesas sit 
atop a superconducting substrate, (a) ^oi is fixed as indi- 
cated, and Bi/G^ii — 7. (b) ^oi is random and Bi/G^ii = 1. 

current density ex Iq and B\\^ac{t) = beyond the skin 
depth (^ 0.15 jam ) inside the BSCCO substrate due 
to the ac Meissner effect, the Ampere-Maxwell bound- 
ary condition forces i^||^ac(^) = just above the BSCCO 
substrate [5]. This corresponds to a perfect magnetic con- 
ductor (PMC) substrate, with the effective image sources 
sketched in Fig. 5 [3, 4]. Thus, for a BSCCO substrate, 
we restrict ^ to < ^ < 90° and replace r]{z') in Eqs. 
(28) and (8) by 

r]-{z') = sgn{z')e[h^ - {z')% (39) 

For cylindrical mesas in the radiation zone, h << a^r, 
and we assume h l/Zc^ for the relevant n. One may ex- 
pand e^^'^-^/R in the generalization of Eq. (29) for small 
z'. From Eq. (39), the electric and magnetic current 
substrate factors are simply obtained from terms linear 
in z\ In the radiation zone, 

S;/l{0) -iknhcosO 0(90'' -0), (40) 

and S^p{6) is obtained from S:/^{6) by kn kmp^ where 
(mp) applies either to the cavity mode, which can be 
either cylindrical or rectangular, as disccussed in Sec. 
IX. For the case of the tail of a broad cavity {mp) mode 
overlapping the nth ac Josephson harmonic, we replace 
S^p{0) with 6'^ ((9), as in Eq. (56) in the Appendix. 

In Figs. 10 and 11, we show our predictions for the 
conducting dipole model with n = 1,2 respectively, of a 
cylindrical mesa atop a superconducting substrate. Note 
that in each case, the radiation vanishes at ^ = 90°, un- 
like any of the radiation from the same modes shown in 
the left panels when the mesa is suspended in vacuum. 
The substrate factor makes a dramatic effect for ~ 90°, 
causing it to vanish at ^ = 90°, in apparent agreement 
with experiment. In Fig. 10(a), we assumed 0oi was 
fixed, and Bi/g'^^i = 7, which maintains the point kink 
found for the combined output pictured in Fig. 6(a) for 
the mesa suspended in vacuum. In Fig. 10(a), the pre- 
dicted radiation vanishes at 6> = 90°, but it is still highly 




FIG. 11: (Color online) Plot of the radiation power at the 
n — 2 second harmonic of the ac Josephson current for a 
cylindrical mesa with k2a = 2* 1.8412/nr on a superconduct- 
ing substrate. 

anisotropic for 9 < 90°, with a point kink that is off 
center (not at ^ = 0°). Experiments so far cannot elim- 
inate this possibility, because scans across the diameter 
of the cylinder have only been made on two perpendicu- 
lar planes, and it is possible that such a point kink would 
have been missed. In Fig. 10(b), Bi/G^i = 1, (j)oi is ran- 
dom, the point kink is removed, and the output power has 
a maximum at ~ 34.5°. These predictions are remark- 
ably similar to the experimental results [26]. In Fig. 11, 
we showed the modification to Fig. 2(b), the second ac 
Josephson harmonic at n = 2, due to the superconduct- 
ing substrate. Since the (01) cavity mode is significantly 
far removed from this ac Josephson mode, this should 
represent the output power a the second harmonic. It 
has a maximum at ^ 42.2°. 



IX. RECTANGULAR MESAS 

Finally, we come back to the more complicated rect- 
angular mesas. A sketch of the electric and magnetic 
surface current densities for the (10) rectangular mesa 
cavity mode is shown in Fig. 12. Although the surface 
electric current density Js is uniform along the edge of 
the mesa, the magnetic surface electric current changes 
its direction by 90° at each of the corners. Hence, a pre- 
cise analytic treatment of 5Jn{x') within the rectangular 
mesa, subject to the correct boundary condition that the 
tangential component H\\ of the magnetic field vanishes 
on each of the mesa edges, is complicated by the presence 
of the corners. However, we shall ignore such complica- 
tions, and focus of the radiation that arises when Js and 
Ms are placed on the mesa edges. 

We now consider a rectangular mesa of width length 
^, and height h in vacuum satisfying ^ > \c > where 
Ac ~ 170 jiia is the c-axis penetration depth. In this case, 
the rectangular mesas allow Josephson vortices parallel 
to the mesa widths, suggesting that the only stable cav- 
ity modes that could be excited are the TM^qq modes, 
which are constant in position along the mesa height and 
length, and oscillate in position with integral multiples of 
half- wavelengths along the mesa widths. In rectangular 
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FIG. 12: (Color online) Sketch of a rectangular mesa with a 
surface electric current density Js (red vertical arrows) and 
a magnetic current density Ms (blue horizontal azimuthal 
arrows). The =b signs refer to the signs of the rectangular 
(10) cavity mode, separated by the yellow line. 

source coordinates {x\y\z')^ 

Js{x') = z'^rij{z')J2[Ux',y')+ga{x',y')m 

<T=± 

Ms„{x') = ^^T]M{z')sm[n{x' - Xn)n/w] 

X E ^iv'M^'^ y') - ^'9a{x', y% (42) 

<T=± 

U{x', y') = w5ix' + aw/2)e[ie/2f - (y')% (43) 

gM, y') = my' + ae/2)e[{w/2r - {x'f], (44) 

where the TM^qq cavity mode energy is degenerate for 
—w/n < Xn < w/n. 

Since standard cavity mode boundary conditions may 
be altered when the corners are treated correctly, we 
treat the output power obtained from the combination 
of the Js and Msn sources in three models. In the 
first model, we treat the combination as coherent, sat- 
isfying Hyix' = ±tL'/2) = 0. This leads to the solu- 
tions Xn = O^w/n for n odd, and Xn = :^w/2n for n 
even. These cases lead to an asymmetry and kinks in the 
angular distribution of the output power. In the other 
two models, the output from the two sources is incoher- 
ent. In Model I, we average the output power P{xn) 
over the two waveforms with the same boundary con- 
ditions, {P{xn))i = l[P{0) P{w/n)] for n odd, and 
{P{xn))i = l[P{w/2n) ^ P{-w/2n)] for n even[3]. This 
model assumes the standard cavity model applies, but 
for mixed sources, we assume that both possible phases 
(such as the signs in Fig. 12) of the electric field are 
allowed with equal probability. This is the rectangu- 
lar cavity analogue of the random 0on niodel for the 
emission from a cylindrical cavity. In Model II, we let 
{P{xn))ii = j^^i^dXnP{xn)n/2w. This model relaxes 
the standard cavity prescription, but preserves the wave- 
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lengths across the mesa. 

Then, the time-averaged power per unit sohd angle in 
the radiation zone is 



dP 



dVt '•/«^°° 1287r2 



UmeanXn{0A)si{e)\ 



+an{9) (Ci + - sin^ 9[Cl cos^ 



+DI sin^ 



El sin 6 cos ( 



(45) 



where i = I, II, an{0) = \EonS^ {0)\y{2ZoJj)^ and 
Xn(<9,(/>) and the C;((9,^), Di^{0,(l)), and Ei{0,(l)) are 
given in the Appendix. Our results for the angular de- 
pendence of the output power from rectangular mesas is 
shown in Figs. 13-15. These figures are for rectangu- 
lar mesas with the length normal to ^ = 0. Although 
the figure "boxes" appear to have the short length along 
^ = 0, that is the direction in which the radiation is 
largest. Since these figures correspond to the (mO) modes 
for m = 1,2, the radiation is primarily along the length of 
the mesa, along or near to ^ = 0. In the figures that fol- 
low, we assume i/w = 20/3, as in some of the mesas [1, 2]. 
Three-dimensional plots of /(^, (j)) cx dP{9^ <p)/dQ in arbi- 
trary units are then obtained. First, we present the data 
for the emission from the primary source, the ac Joseph- 
son current in the form of the surface electric current den- 
sity Js- In Fig. 13, we present I{0, (j)) for the n = 1 and 
n = 2 ac Josephson radiation, respectively for the mesa 
suspended in vacuum. The corresponding predictions for 
the rectangular cavity (10) and (20) modes at the same 
frequencies as the n = 1,2 ac Josephson radiation are 
shown for Model I in Fig. 14 for the cavity suspended 
in vacuum. We note that the output for the rectangular 
cavity (10) mode pictured in Fig. 14(a) is very similar to 
that for the cylindrical cavity (11) pictured in Fig. 4(a) 
when both mesas are suspended in vacuum. We note that 
the angular dependence of the second harmonic emitted 
from rectangular cavities should be very interesting to 
measure, especially as it contains a mixture of the radia- 
tion predicted in Figs. 13(b) and 14(b), and is distinctly 
different in form from any of the cylindrical cavity modes 
calculated without a substrate. 

In Fig. 15, we show the combined output from the 
primary ac Josephson current and excited cavity modes 
for a rectangular mesa atop a superconducting substrate. 
In Fig. 15(a), the fundamental n = 1 mode is locked onto 
the rectangular cavity (10) mode, and the figure assumes 
a(0) = 0.2 using Model I for the cavity. This figure is 
very similar to that observed experimentally [24]. In Fig. 
15(b), the predicted output power is shown for the second 
harmonic, again using Model I with a(0) = 0.2 for a mesa 
atop a superconducting substrate. 

As for cylindrical mesas, the output from rectangular 
mesas is linearly polarized. Because the angular distribu- 
tion of the output intensity at the fundamental frequency 




FIG. 13: (Color online) Three-dimensional plots of the radi- 
ation intensity in arbitrary units for rectangular mesas with 
l/w = 20/3 from the uniform ac Josephson current alone, 
when the mesa is suspended in vacuum, (a) At the funda- 
mental n = 1 frequency with kiw — Tv/rir. (b) At the second 
harmonic n — 2 with k2W — 2ti /ur. 



0=0 




FIG. 14: (Color online) Three-dimensional plots of the radi- 
ation intensity in arbitrary units for rectangular mesas with 
i/w = 20/3 from the cavity source alone, when the mesa is 
suspended in vacuum, (a) The fundamental (10) cavity mode 
with kiow = Tv/rir. (b) The second harmonic (20) cavity 
mode with k2ow — 27v/nr. 



is very similar in both rectangular and cylindrical mesas, 
only minor differences in the angular distribution of the 
tilt angles for the two types of mesas are expected. 

We now consider the case when the combination of the 
radiation sources is coherent. In this case, the output 




FIG. 15: (Color online) Plots of the intensity in arbitrary 
units for rectangular mesas with i/w = 20/3 of the combined 
radiation from Model I with a{0) — 0.2 from the uniform ac 
Josephson current and the rectangular cavity modes when the 
mesa sits atop a superconducting substrate, (a) The funda- 
mental n = 1 ac Josephson mode and the cavity (10) mode 
with kw = kiow = Tv/rir. (b) The second harmonic with 
k2W = k2ow = 27v/nr. 
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FIG. 16: (Color online) Plots of the intensity in arbitrary 
units for rectangular mesas with ilw — 20/3 of the com- 
bined radiation when the mesa sits atop a superconducting 
substrate, (a) The fundamental n — \ ac Josephson mode 
and the cavity (10) mode with kw — kiow — Tr/rir, xi — 
and ri — 0.05 (b) The second harmonic with k2W — k2ow — 
27r/nr, X2 — w/4: and r2 — 0.6. 

power is given by 
dP 

— (X |-sin(9xn+rn(M^sin(^-M^cos^)P 

+r^|M;^cos(/) + M^sin(/)p, (46) 

where = EquS^ {0)/[2iiioJjancS^ {0)] is real. For 
n = 1, Mf and , given in the Appendix, when evalu- 
ated at either xi = or w are real, and a coherent combi- 
nation is strongly asymmetric. For n = 2, however, both 
Mf and are pure imaginary for X2 = Hence, 
the resulting coherent combination of the output power 
is a sum of the two output powers from the two sources. 
This is pictured in Fig. 16. In Fig. 16(a), we show the 
combined coherent output when ri = 0.05, computed 
when xi = 0. This figure shows a strong intensity asym- 
metry even for this small ratio of the outputs from the 
two sources, with the pattern exhibiting Ci symmetry. 
In Fig. 16(b), the combined output for n = 2 is shown 
for the case r2 = 0.6, computed with X2 = w/A. In this 
case, the pattern has C2v symmetry. 



X. DISCUSSION AND SUMMARY 

In this paper, we have presented the first correct pre- 
scription of how to incorporate the Ampere boundary 
condition into the superconducting stack of Josephson 
junctions under the application of a dc voltage V, which 
itself generates THz radiation. We have used both of 
the Love equivalence principles. By the Love equivalence 
principle for a magnetic conductor, the radiation from 
the magnetic field inside the mesa generated by the ac 
Josephson current is effectively set equal to zero inside 
the cavity, and replaced by the equivalent surface elec- 
tric current J^, which is the primary radiation source. 
We remark that setting H = inside the mesa also sets 
the boundary condition for the cylindrical cavity. With 
this boundary condition, the cylindrical cavity modes are 
immediately found, and when the ac Josephson current 



effectively placed on the edge as Js radiates at a fre- 
quency that matches one of the cavity modes, the inho- 
mogeneous part of the ac Josephson current has a mode 
which resonates with a cavity mode of the same (or per- 
haps similar, in an extension of our model) spatial form, 
and the amplitude for that cavity mode then grows lin- 
early in time, until it saturates, and then radiates in con- 
junction with the primary ac Josephson radiation at that 
frequency. 

The form of the cavity radiation is obtained from 
Love's electric conductor equivalence principle, in which 
the equivalent of the cavity electric field is placed on the 
surface as a magnetic surface current density Ms^ which 
radiates in conjunction with the ac Josephson radiation 
in Js, the inhomogeneous part of which excited the cav- 
ity resonance. From the experiments on three cylindrical 
mesas, we conclude that the lowest energy cavity (11) 
mode has been excited. In addition, the harmonic radia- 
tion at twice and three times the fundamental frequency 
cannot be explained from the cavity alone, as the mis- 
match in frequencies is 4% for the second harmonic, and 
about 3% for the third harmonic. This harmonic radi- 
ation therefore arises almost entirely from the primary 
radiation source, the ac Josephson radiation acting as a 
surface electric current density. In addition, the fact that 
the superconductor acts both as an electric conductor 
and as a magnetic conductor provides a mechanism to 
understand the role of a superconducting BSCCO sub- 
strate. Such substrates cause the emitted radiation to 
vanish along the direction parallel to it, greatly reducing 
the output power of the radiation. 

We remark that while the main function of the inho- 
mogeneous ac Josephson current is that it couples to a 
cavity mode, provided that the conditions are properly 
met, if the inhomogeneities are sufficiently strong, they 
can significantly alter the azimuthal and axial anisotropy 
of the combined radiation at the fundamental frequency, 
and especially that of the higher harmonics. In addition, 
it remains an open question as to whether the angular 
constants (j)on are fixed or random during the time of the 
measurement. To the extent that the cylindrical cavity 
itself is perfectly homogeneous, one might expect no pre- 
ferred angle for each mode of inhomogeneity. On the 
other hand, some feature involved in the experimental 
situation might lead to a fixed inhomogeneity direction. 
Only experiments measuring the angular distribution of 
the radiation and its polarization can distinguish these 
two situations. Although we did not show any figures to 
illustrate this point, if one were to assume the phase of 
the cavity mode for a fixed 0on could have either sign 
with equal probability, then the situation would be anal- 
ogous to our Model I for a rectangular mesa. For a cylin- 
drical cavity, this model would preserve the point kink of 
the output radiation and restore C2v point group symme- 
try, but the position of the kink would be at ^ = 0°, and 
the value of the output power would not vanish there. 
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For rectangular mesas, setting n x H = on each of 
the edges leads to singularities in the analytic properties 
of the mathematical difficulties at the corners, greatly 
complicating the situation. Hence, most previous work- 
ers studied infinite strips, where the corners could be ne- 
glected. For the case of finite length rectangular mesas, 
the question of whether the combined radiation is skewed 
to one side of the mesa or the other is also relevant. If 
the mesa were symmetric, as in a cylinder with no pre- 
ferred angle, one would average over the two configura- 
tions, resulting in a symmetric radiation pattern, and 
an incoherent combination of the radiation from the two 
sources, each of which is separately coherent. That was 
the purpose of our study of Models I and II. 

We remark that we have treated the cylindrical mesa 
essentially as an infinitely thin disk, except in calculating 
the superconducting substrate factor, completely neglect- 
ing any spatial variation normal to the layers. While this 
is not a good approximation for treating heating inhomo- 
geneities, it suffices for impurity inhomogeneities and the 
important lateral inhomogeneities in the ac Josephson 
current, which we assume to be the critical part of the 
microscopic mechanism that excites the cavity modes. 
Hence, exotic features such as inhomogeneous kinks along 
the z axis are omitted entirely, as they are completely 
unnecessary for the excitation of the cavity modes, but 
instead complicate the situation unnecessarily [31]. 

We emphasize that cylindrical mesas are much simpler 
to understand theoretically than are rectangular mesas. 
There are several aspects to this simplicity. Most impor- 
tant, the von Neumann boundary condition is trivial in 
this geometry, whereas for a rectangle, it is practically in- 
tractable at the corners, except by numerical techniques. 
In addition, the fact that none of the cylindrical cav- 
ity modes are harmonics of one another is an extremely 
important point, which is why we set out to study cylin- 
drical mesas in the first place. Since rectangular cavities 
have higher energy modes that are also harmonics of a 
lower mode, such as the (10) or the (01) mode, the obser- 
vation of harmonics in the output of rectangular mesas 
did not allow for a precise determination of the primary 
radiation source, which has now been clearly identified 
as the ac Josephson current, and not the cavity mode 
radiation. The primary role of the cavity mode is to lock 
an ac Josephson mode (usually, if not always, the fun- 
damental) onto a well-defined cavity mode, fixing coj to 
(jOmopo • Then, the laterally inhomogeneous ac Josephson 
current excites the (moPo) cavity mode, allowing it also 
to radiate. From fits to experimental data on both rect- 
angular and cylindrical cavities, it appears that at least 
half of the intensity at the angles ^max ~ 30 — 40° of max- 
imum output arises from the primary source. This pri- 
mary radiation source radiates at all of the ac Josephson 
harmonic frequencies, and to date, the second harmonic 
has been observed in emission from all three cylindrical 
mesas under study, and the third harmonic was visible in 



two of the three mesas studied. The amplitudes relative 
to that of the fundamental of the higher harmonics are 
comparable to those in rectangular mesas, for which am- 
plification by the excitation of higher cavity modes could 
occur. 

The cavity modes of a cylindrical mesa are similar in 
form to those of a drum, and one or more of them can 
be amplified when a drummer strikes a particular spot 
on the drum surface. The striking surface region of the 
drum does not have to have the precise shape of that of 
the main cylinder (drum) mode excited. By analogy, a 
modification of our theory of the cylindrical cavity mode 
amplification could occur without implementation of the 
von Neumann boundary condition to provide a precise 
matching in spatial form of the cavity mode with that 
of the inhomogeneities. In a rectangular mesa, a subset 
of the cavity mode frequencies are harmonics of one an- 
other, and can be excited as by a player of a stringed 
instrument in lightly touching a finger at the midpoint 
or quarter point at either end of the string, for example, 
in order to make it sound one or two octaves higher than 
the fundamental, respectively. 

Finally, we note that this is the first treatment of the 
dramatic effect of superconducting substrates upon the 
output power of radiating BSCCO mesas. Using the 
same Love principle that was the basis for the correct 
implementation of Ampere's law, the magnetic equiv- 
alence principle, we conclude that the superconducting 
substrates cause a drastic reduction in the output power 
of the radiation, especially for the power emitted near 
to ^ = 90°. Note that in our predictions for both the 
primary ac Josephson current and the secondary cavity 
radiation sources, the output power at the fundamental 
frequency is predicted to be at least comparable to that 
at its maximum output. The fact that in both rectan- 
gular and cylindrical mesas, the experimental output at 
= 90° is consistent with zero is a very strong indication 
that our magnetic conductor model of superconducting 
substrates is correct. Our model is also consistent with 
the experiments of Barbara et a/., in which the output 
of a Josephson junction array with the currents parallel 
to the superconducting substrate was enhanced by the 
substrate (or ground plane) prior to its entry into the 
waveguide [3 2-34]. We reiterate that removal of the su- 
perconducting substrate could enhance the output by at 
least two orders of magnitude, and by replacing the sub- 
strate with a perfect electric conductor such as Cu, one 
could further enhance the output by a factor of four. This 
could allow for output as high as 5mW, which would be 
more than sufficient for many practical applications. 

In summary, we have identified the primary micro- 
scopic source of the coherent radiation as the ac Joseph- 
son current which locks onto a cavity mode, exciting it, 
and the two modes radiate together at the fundamental 
ac Josephson frequency. It suffices to treat all of the junc- 
tions as acting in unison. For radiation at the fundamen- 
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tal ac Josephson frequency locked onto the cylindrical 
cavity (11) mode, we predict a combined output radia- 
tion pattern with linear electric field polarization that is 
usually along 0^ but not always. The second harmonic 
in cylindrical mesas arises almost exclusively from the 
ac Josephson current, and is also linearly polarized, the 
output from rectangular mesas should also be linearly 
polarized. We reiterate that removal of the supercon- 
ducting substrate, or better yet, replacement of it by a 
perfect electric conductor such as Au or Cu, could lead to 
an enhancement of the output power of the mesas from 
the highest observed power of 5/i W up to 5 mW, suitable 
for many applications. We therefore name this device a 
Josephson STAR-emitter, for stimulated terahertz am- 
plified radiation emitter. 
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APPENDIX 
A. Electric field from 5Jn{x') 

We now calculate the electric field EAsj{x,t) arising 
from the 6Jn{x) contribution to the electric surface cur- 
rent density. From Eq. (10), we write, 

Asj{x,t) = d^x'z'n{z')5{p' -af- 

n=l 

oo 

^ Yl C^mpJm{k'mpp')cOs[m{4>' -(f>On)]. 

p=l,m=0 

(47) 

In the radiation zone, we have 

CO 

Asj{x,t) -Osinepov Yl C^r^^i-^r 

n,p=l;m=0 
^in{kjr—(jjjt) 

X Jm {k'mp(^) Jm {nke ) 

xcos[m{(^-(^^n)]Si{e). (48) 



We then find 

CO 

n,p=l;m=0 
^in{kjr-ujji) 

X Jm {k'^po) Jm {nke ) 

xcos[m(0-0on)]^;!(^). (49) 
Thus, combining this with Eaj{x, t), we have 

n=l 

CO 

xnojjS:/,{0) Y(-iT^^EmnJm{nke) 

m=0 

X cos[m(0 - 0on)], (50) 

oo 

Emn = ^n^m,0 + ^C'^"^^m(Cp«). (51) 
p=l 

The Ejnn are real. Thus, there are two effects, both 
of them we assume to be small, of including the full 
anisotropy of the spatial dependence of the ac Josephson 
current into the surface electric current density source 
Js. The first is a renormalization of the uniform co- 
efficients to 

CO 

Jn = Jn+T.^^Po{Kp^). (52) 
p=l 

and the second is to introduce azimuthal anisotropy 
and additional axial anisotropy in the additional terms 
Jm{nke) cos[m{(f) — 0on)]- We note that the renormal- 
ization in Eq. (52) is a bit spurious, as the term would 
vanish if the integral were over the full volume of the 
cylinder, as noted following Eq. (20), at least if the phase 
factor e~^^'^ were absent, but we have placed the cur- 
rent only on the edge, modifying the integral. The second 
change is more interesting. It means that output power 
which is azimuthally anisotropic can arise from inhomo- 
geneities in the ac Josephson current itself, and does not 
necessarily depend upon the presence of the cavity. How- 
ever, we generally assume the Cmp <^ Jn^^^ that the pri- 
mary effect of these spatial inhomogeneities 5J^{x') is to 
excite a particular cavity mode, as discussed in Sec. IV. 
Nevertheless, as discussed in Sec. VI, the case of n = 2 
is of particular interest, as Eq. (50) demonstrates that 
a ^-dependence to the output power can arise from the 
inhomogeneous ac Josephson current, although for every 
n, the output power from such terms is proportional to 
sin^ ^, vanishing at ^ = 0°. As we shall see, if the n = 1 
ac Josephson mode locks onto the (11) cavity mode, the 
nearest cavity mode to the n = 2 ac Josephson harmonic 
is the (01) mode, which should give rise to a very weak 
output, at best, but one that would be ^-independent, as 
well. Another point of some interest is that the mth con- 
tribution to Ea is proportional to (—1)^+-*^. This could 
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be important when the radiation originates both from the 
inhomogeneous part of the ac Josephson Js and cavity 
Ms sources, as discussed in Sec. VI . 



B. Resonant and ofF-resonant cavity radiation 

The F vector potential resulting from the resonant cav- 
ity mode is then 



Jm{k'mpa) — 0'cos[m(0' -0on)] 



xe 



-ikmpOL sin 9 cos((f)—(f)') 



where cori 



nix) J and /c^, 



(53) 



Jc outside the 



mesa, and S^p{0) = 1 when the mesa is suspended in 

vacuum. After setting (j)' = ^cos^sin(0 — 0') + (/)cos(0 — 
(/)'), and evaluating the integral over 0', we obtain 



F(x,t) 



r/a»l 



-(-if 



krr 



47rr 



^cos6>sin[m((/) - (t>om)]J^{k^p) 

-<^COs[m(0 - (t)On)]J^{k^mp) 



where 



(54) 



^mp = kmpasinO, and J^{z) = [J^+iW ± Jm-i{z)]/2. 
The same procedure can be followed for the F vector 
potential arising from the off-resonant co = nujj harmonic 
electric current source exciting the tail of an (mp) cavity 
mode, where {mp) ^ {mopo), and Eqs. (54) and (55) are 
again obtained, with the replacements that ujmp tiujj, 
kmp nkj = nujj/c in Eq. (54), and Gmp -iGmp, 
where 



mp 



Jm{k'a)S^{0), (56) 



where the substrate factor {0) = 1 when the mesa 
is suspended in vacuum, and otherwise differs only from 
S^p{0) in that kmp nkj^ as noted in Sec. VIII. 



C. Vector potentials and surface currents for 
rectangular mesas 

For the rectangular mesas studied in Sec. IX, A(x,t) 
and F{x,t) in spherical coordinates are then given in the 



radiation zone by 



A{x,t) 



Snr 



^a„e-(*'-'— -*)5;((^)Xn,(57) 



n=l 



Xn{0,(t)) = cosXn— hcosFn^; — , (58) 



F{x,t) 



Mi 



My 



^^One'"('=-"^-'^-'*)5f (6») 



x{xM^+yMy), 



(59) 



-siny„ Yl (Te™"/-?^^fe^, (60) 



^e-^^sin(^ + ^)(61) 



where = (A:^it'/2) sin^cos(/), = (A:^-^/2) sin^sin(/), 
knW = nTr/rir, ^n,cr = n7r/2 -h crXn, V = w£h, 
S^{0) = S^{0) = 1 for no substrate, = r sin 6> cos -h 
6 cos cos — sin and ^ = r sin ^ sin -h ^ cos sin 
(/)COS(/), and Xn appears in Eq. (42). 

The quantities present in the expression for the output 
power intensity from the combined electric and magnetic 
surface current density sources follow. For the portion of 
the output power resulting from the uniform portion of 
Js: Eqs. (57) and (58) are sufficient. The part of the 
combined output arising from Ms is more complicated, 
and the combination can be either coherent or incoherent, 
as for the output from cylindrical mesas. When the com- 
bination is coherent, we assume that Hy{x' = ±.w/2) = 
is maintained along both lengths of the mesa, but the out- 
put power P{xn) is evaluated when Xn is either ^ oi w/n 
for n odd, or when Xn = -^w/2n for n even. When the 
output from the two sources is incoherent with respect 
to one another, then we average P{xn) in two models. 
In Model I, we also assume that Hy{x' = ±.w/2) = 
is maintained along both lengths of the mesa, so that 
{P{xn))i = l[P{0) + P{w/n)] for n odd, and {P{xn))i = 
^[P{w/2n) + P{—w/2n)] for n even. In Model II, we re- 
lax the boundary condition upon Hy , and average over all 
Xn values that preserve the wave vector within the mesa, 
{P{xn))ii = {n/2w) J^l'j^P{xn)dXn. In both models, 
the average output power is characterized by three func- 
tions Cl{e,(t)), Di^{0,(t)), and ^;(<9,0) for i = I and II. 
For Model I, d = Al, = B^, and = 2AnBn, 
where for n either odd or even. 



An = sinYn ^ 



sinl^ . fUTT 



■ sm 



For Model 11, 



(62) 
(63) 

(64) 
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D. 



y2 (l-(-ircos(2X0), (65) 

(66) 



^11 ^ p Sin^^n ^ crsin^(X^,^) 



cr=± 



where X^, Y^, and X^^a given following Eq. (61). 
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